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ABSTRACT 

We present a detailed statistical analysis of the alignment of polarizations of radio 
sources at high redshift. We use the JVAS/CLASS 8.4-GHz surveys for our study. This 
study is motivated by the puzzling signal of alignment of polarizations from distant 
quasars at optical frequencies. We explore several different cuts on the polarization 
flux for our analysis. We find that the entire data shows a very significant signal 
of alignment on very large distance scales of order 500 Mpc. The alignment starts 
to decay only at much larger distances of order Gpc. If we only consider data with 
polarization flux greater than 1 mJy, we find alignment at distance scales less than 150 
Mpc. We also find that data with polarization flux less than 0.5 mJy does not show 
significant alignment. Similar results are seen for data with degree of polarization less 
than 0.01, although here a mild signal of alignment is observed for a narrow range of 
angular separations. We argue that the signal cannot be explained in terms of bias 
due to systematic errors in removal of instrumental polarization. We also find that the 
degree of polarization shows a strong negative correlation with the total flux. The data 
appears to fall into two classes, one of which shows such a correlation. The remaining 
set, which has total flux greater than 100 mJy and degree of polarization lying between 
0.01 and 0.1, shows a more random behaviour. The latter set is also found to show no 
alignment whereas the first set shows a very strong polarization alignment. 
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1 INTRODUCTION 

The Big Bang model assumes that the Universe is homoge- 
neous and isotropic on large distance scales. However there 
currently exist several observations which appear to violate 
this basic assumption. In particular the radio polarizations 
from radio galaxies a ppear to show a large scale dipole pat- 
tern across the sky (IBirchl Il982l ; iKendall and Young! 1 1984 
Ijain and R alston 1 9991 ). The possibility that the signal ob- 
ser ved by IBirchl (I1982T) might a rise due to bias was raised 
in (jPhinnev and Webster 19831). Th e signal was dismissed 
by iBietenholz and Kronberj fl984 ) who found that it is 
not pr esent in a larger data set. However I Jain and Ralston] 
l|l999l ) argued that the signal is present if we consider all 
the radio sources for which the relevant information, i.e. 
the polarization position angle and the galaxy orientation 
angle, is available in the li terature. It is also inte resting 
that the dipole axis found in I Jain and Ralston! l| 19991 ) aligns 
closely with the Cosmic Microwave Background Radiation 
(CMBR) dipole axis. Furthermore a recent study finds a 
dipole anisotropy in the brightness of radio sources which is 
much larger than what is pred icted by the Doppler effect due 
to local motion (|Singalll201ll ). This also indicates a dipole 



axis which may be of cosmological origin and well aligned 
with the CMBR dipole. 

The optical polarizations from quas ars also show align- 
ment over cosmologically large distances (lHut scmckcr s~et alj 
Il998l : iHutsemekers and L amv 2000l; I Jain et aL I12003J). The 
distan ce scale o f these correlations is of the order of 
Gpc |jain et all 120031 ). Furthermore the Cosmic Mi- 
crowave Back ground Radiation (CMBR) d a ta shows sev- 
eral fe atures (jd c Olivcira-Costa et al. 2004: Eriksenljtjil] 
2004 iRalston and Jainl 12004 lLand and Magueiiol 120051 : 



Kim and NaselskvT 1201 
not consistent with 



Samal et aU 200: 



20091 ) which 
We point 



are not consistent witn Big Bang cosmology, 
out that the WMAP science team has argued that some 
of the claimed anomalies in CMBR data set may arise 
due to a posteri ori choice of statist ics to test for a 
particular effect (|Bennett et al.l 1201 ll ). In other words 
one notices a particular odd feature in the data and 
then devises a statistic to test its significance. Such a 
procedure is likely to overestimate the significance of 
the detected anomaly. In view of this it is extremely 
interesting th at alignment axis of CMB R quadrupole 
and octopole Jde Oliveira-Costa et al.l 120041 ) . radio dipole 
axis (|jain and Ralston! 19991 ) and the two point correla- 
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tions in the optical polarizat i ons dHuts cmckcrs et al.|[l998l ; 
IHutsemekers and Lamy||2000l ; |jain et al.ll2003l) all align very 
closely with the CMBR dipole axis ( Ralston and Jainll2004h 
and point roughly in the direction of the Virgo cluster. Fur- 
thermore there has also bee n claims of violation of is otropy 
in cluster pecul iar velocities dKashlinskv et~a l, 2009b aj) an d 
galaxy surveys (|ltoh et al.ll2010l ). Remarkably the cluster pe- 
culiar velocities also indicate a direction close to the CMB 
dipole if we include the highest redshift data. The spiral 
galaxies a lso appear to show an interesting signal of parity 
violation (|Longoll201ll ). There exist many attempts to the- 
oretically explain these observations, most of which assume 
violation of the cosmological principle. An interesting pos- 
sibility, which is completely consistent with the Inflationary 
Big Bang model, is that the universe was inhomogeneous 
and anisotropic at very early stage, before the epoch of infla- 
tion. It evolves into a homogeneous and isotropic de-Sitter 
space-time during inflation. It has been shown that, for a 
wide class of models, there exists parameter range such that 
the modes generated during this early phase can re-enter 
the horizon much before the current era and h ence can af- 
fect present observations (|Aluri and Jain|[201ll ). These can, 
in principle, generate the observed anisotropies. 

In (| Jackson et al.l [20071 ) . the authors compiled db CcI.1jcL _ 
logue of radio polarizations from distant radio galaxies. The 
catalogue contained a total of 12743 sources. Motivated by 
the observed linear polarization alignment in quasar d ata 
l|Hutsemekers et al.lll998l; IHutsemekers and Lamvll2000h at 



visible wavelengths, iJoshi et al. |2007l ) studied the possi- 
bility of a similar effect at radio frequencies using this 
catalogue. No significant effect was detected. The align- 
ment in optical polarizations is seen over cosmologically 
large distances of order Gpc (IHutsemekers et al.l 19981 ; 



jp c , 

IHutsemekers and Lamvl 120001 ; Ijain et al" 20031 ) . and the 



phenome na is very puzzling. Th e re are several possible 
models (IHutsemekers et al.l 19981; IHutsemekers and Lamv 



200d ; |jain et al.ll2002l; IHutsemekers et alj|2005l; IPavez et al 



2008; iPiotrovich et all 120081 ; lUrban and Zhitnitskvl l20ll 
Ciarcellutill2012D looking for the interpretation of this align- 



ment. One possibility is that this alignment arises due to 
large scale correlatio ns in the intergalactic magnetic field 
(|Agarwal et al.|[2008l ). The inter galactic magnetic field may 
be seeded in the early Universe dSubramanian and Seshadril 
l2003l ; ISeshadri and Subramanianll2005l , 120091 ) and presence 
of large scale correlations in this field are not neces- 
sarily in conflict wi th the Big Bang cosmological model 
|Agarwal et al.|[2008l ). We still require a mechanism for how 
such a magnetic field generates large scale correlations in the 
optical polarizations. One possibility is that this is caused 



dJain et al.ll2002 


Hutsemekers et al.ll2005l;|Pavez et al.ll2008l; 


IPiotrovich et al 


20081; Agarwal et al.l 20081). This effect is 



freq uency dependent and is much smaller at radio frequen- 
cies (|jain et al.ll2002D . Hence it may be consistent with ab- 
sence or reduced alignment effect in radio polarizations. 
There have also been other proposals to expl ain this effect 
(lUrban and Zhitnitskvll20Tll ; ICiarcellutill2012t ). 

In the present paper we analyze the radio data in or- 
der to study possible alignment of radio polarizations. A 
study of radio pola rizations to check for alignment was ear- 
lier conducted in (|joshi et al.ll2007l ). Here we extend this 
work by considering the fact that polarization angles de- 



pend on the coordinate system used. In order to properly 
analyze the presence or absence of large scale alig nment one 
needs to define a coordinate invariant statistics (|jain et al.l 
l2003f ). The basic point is that, in order to compare two po- 
larization angles on the surface of the celestial sphere, one 
needs to parallel transport one of them to the position of 
the second along the great circle joining the two points. The 
contribution due to parallel transport may be negligible in 
most cases, if the two points are separated by a small dis- 
tance. However it becomes very important if we are testing 
align ment over large distances. Furthermore in (|joshi et al.l 
120071 ). the authors restricted their analysis to data which 
has polarization flux density greater than 1 mjy. Here we 
also analyze data with lower polarization flux density. 

The data at low polarization flux density might be 
more contaminated with noise. However this will not af- 
fect our alignment results unless there is significant bias in 
the data. A potential source of bias is the e rror in the re- 
mova l of residual instr umental polarization (|jackson et al.l 
120071 ; IJoshi et al.ll2007l ). This can lead to large scale cor- 
relations in polarizations even when none are present. It 
is clear that this effect will dominate for sources which 
have low degree of polarization and/or low polarization 
flux. Hence one can evaluate its contribution by focussing 
on data with lo w polarization. Another source of error is 
the p ositive bias (|Simmons and Stewart|[l984l ; I Jackson et al.l 
120071) that arises in the degree of polarization. This arises 
since the degree of polarization depends on the squares of 
Q and U and will always acquire a positive value. How- 
ever this cannot affect the alignment statistics which are 
based only on the linear polarization angle. Furthermore 
we point out that there is some motivation for looking at 
alig nment at lower p o lariza t ions since some theor etical mod- 
els (|Das et all 120051 . |200S| ; lAgarwal et al. I |200S| ) predict a 
smaller polarization for radio frequencies. Hence here we 
might expect alignment only at smaller polarization. We 
also point out that the entire field of Cosmic Microwave 
Background Radiation polarization is based on polarization 
flux which is of order 10~ 6 in compari son to the total inten- 
sity (|see. for example. We inberg 200^). This is much smaller 
than the degree of polarization of the radio data we consider 
here. Hence we do not see any reason to discard data with 
smaller polarization flux. In any case one should be cautious 
in interpreting our results for low polarizations. 

We also point out that errors in polarization position 
angle (PA) calibration can also reduce the true significance 
of alignment that might be present in data. This will lead 
to a systematic error in the observed PA. This can mask an 
alignment effect present in data si nce this systemat ic error 
may be different in different runs (|joshi et al.ll2007l) . 



2 DATA SELECTION 

In our study we use the data av ailable in the catalogue 
produced by Ijackson et al. I l|2007f ). It contains a total of 
12743 radio sources and lists the angular positions and the 
Stokes /, Q and U parameters. Since the redshift of most 
of these sources is unknown, we assume that these sources 
are roughly at the same redshift equal to unity. The input 
observable for the alignment study is the polarization PA. 
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The calibration methods and the catalogue pro duction has 
been discussed in detail in I Jackson et al.l |2007l ). 



3 STATISTICAL PROCEDURE 

All astronomical observations are made on the hypotheti- 
cal celestial sphere and any directional measurement on this 
sphere correspond to a particular coordinate system. For ex- 
ample, the polarization angles are measured in a local frame, 
formed by two unit vectors <j) an d 0. These unit vectors de- 
pends on the coordinate system used, i.e. they depend on 
which direction we choose as our North Pole, and hence one 
cannot directly compare vectors at different positions on the 
celestial sphere. The proper procedure to compare such vec- 
tors is to transport one of them to the position of the second 
along the geodesic joining t he two positions. A detailed pro- 
cedure has been discussed in ljain et all (12003! ) and we follow 
this procedure for comparing polarization angles of two dif- 
ferent sources. 

We define a statistic Sp to quantif y the alignment 
( Hutscmc kers et al. Hl998l ; I Jain et al.ll2003h . Let us suppose 
we want to calculate the alignment of the polarization at a 
site k with its n v nearest neighbours. We first need to par- 
allel transport the polarization angles from each site to the 
site k and then we compare all polarization angles with po- 
larization at the k th site. We define a measure of dispersion, 
dk, at site k as 



d k = — ^ cos[2(V>; + Ai^k) 
n.~. * — ' 



(1) 



Here ?/>; are the polarization angles of the nearest neighbours 
of the k th site and the factor Ai_>fc arises due to parallel 
transport from i — > k. The sum in Eq. [1] includes the site k 
also. We next maximize dk as a function of ipk- The resulting 
value of ipk is interpreted as the mean polarization angle at 
the site k and the corresponding maxima of the function 
in Eq. [1] gi ves an estimate of dk- The statistic may n ow be 
defined as (jHutsemekers et al.|[l99Sl : Ijain et al.ll2003h . 



1 n s 

Sp = — / dk\ 



(2) 



where n s is the total number of samples in the data. A large 
value of Sp indicates a strong alignment between polariza- 
tion vectors. 

In this paper we use the statistic defined in Eq. [2] to 
determine whether t he radio polarizat ion data shows signif- 
icant alignment. In (|joshi et al.l 120071 ) the authors imposed 
a cut on polarization flux density to include only sources 
with polarization flux greater than 1 mjy. As we discuss in 
the Introduction we study alignment both with and without 
this cut. The data set contains a total of 4400 sources with 
polarization flux > 1 mjy and 7452 sources with polariza- 
tion flux >0.5mJy. The number of sources with polarization 
flux lying between 0.5 mjy and 1 mjy is 3051. The total 
number of sources without any cut is 12743. We note that 
the sources lie dominantly in the Northern hemisphere. Fur- 
thermore there are very few sources along the galactic plane. 

The significance of alignment in the data set is com- 
puted as follows. We first compute the statistics defined in 
Eq.[2]for a given number of nearest neighbours of any source. 




400 600 
Distance (Mpc) 



1000 



Figure 1. The mean distance (Mpc) among sources as a func- 
tion of the number of nearest neighbours (n v ) for the full data 
set (solid line), set with polarization flux greater than 0.5 (short 
dashed line) and with polarization flux greater than 1 (long 
dashed line). 



This statistic is compared with the result for a large number 
of random samples, which are generated by shuffling all the 
PAs among different sources. In our simulations we use a 
total of 1000 random samples for a given number of nearest 
neighbours. The probability that the alignment seen in data 
might arise as a random fluctuation is determined by the 
number of random samples which shows a larger value of 
statistic in comparison to the real data. If we find that none 
of the random samples shows a larger value of the statistic 
then the significance is computed by determining the mean 
and standard deviation of the random set and assuming that 
the distribution of random samples is approximately Gaus- 
sian. 

The number of nearest neighbours n v of any source are 
computed by assuming that all the sources are located at 
the same redshift of 1.0. The redshift information of these 
sources is not currently available and hence we make this 
assumption. Essentially here we are ignoring the third di- 
mension and determining the nearest neighbours only on 
the basis of the angular separations. It is likely that in many 
cases this will lead to a wrong assignment of the set of near- 
est neighbours of a source. However this also cannot gener- 
ate alignment in a data set if none is present. If the data 
does show alignment at any distance scale then it will affect 
the detailed numerical results. For example, let us assume 
that the radio polarizations are aligned over a small distance 
scale of a few Mpc, but show no alignment over larger dis- 
tances. Our nearest neighbour assignment may include some 
sources which are in fact much further away. It is clear that 
these sources which are mis-identified as nearest neighbours 
will only add noise to the signal and reduce the significance 
of alignment. In Fig. 1, we show the relationship between 
the number of nearest neighbours and the mean comoving 
distance from a source within which these nearest neigh- 
bours reside. Here the mean is taken over the entire sample. 
We have assumed that all sources are located at a redshift 
of 1. Furthermore we assume the standard Lambda Cold 
Dark Matter model for computin g the comoving distance 
|see. for example. Weinberell2008l ). 

The distribution of statistics for random sample is 
shown in Fig. 2 for the number of nearest neighbours n v = 10 
and 50. Here we use 1000 random samples generated by ran- 
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Figure 2. The statistic So histogram for the shuffled PA,s 
for 10 and 50 nearest neighbours. The plot shows distribution 
of 1000 random samples generated by randomly shuffling data 
for radio sources with polarization flux lying between 0.5 and 
1.0 mjy. For the case of 50 nearest neighbours, we also show a 
Gaussian fit to the distribution. The verticle line in each plot 
shows the corresponding value of So for the actual data. 

domly shuffling data for radio sources with polarization flux 
lying between 0.5 and 1.0 mjy. For the case of 50 nearest 
neighbours we also show the Gaussian fit to the distribu- 
tion. We find that it provides a good fit to the data. Similar 
distributions are seen for all the values of n„and the cuts 
on polarization flux used in this paper. In all these cases we 
find that the Gaussian provides a good fit to the data. 



4 RESULTS 

In Fig. 3 we show the results for polarization flux density 
greater than one. The upper graph shows the statistic Sd as 
a function of the number of nearest neighbours. The dots de- 
note the values of the statistic for the data sample. We also 
show the results obtained from the random samples gener- 
ated by shuffling polarizations among different sources. We 
show the mean and standard deviation of the random sam- 
ples. The significance of alignment of this data set is shown 
in the lower graph. Here we show the significance in terms 
of the sigma values by which the data differs from random 
simulations. The corresponding results for polarization flux 
density lying between 0.5 mjy and 1 mjy are shown in Fig. 
4. 

For polarization flux greater than unity we do not find 



a significant alignment for number of nearest neighbours 
n v > 14. This is in agreement with the results obtained 
in (jjoshi et al.ll2007ri. However we find significant alignment 
for n v ^ 14. In i Joshi et alj|2007h this region was never ex- 
plored. The significance of alignment is found to be more 
than 3 sigmas for n v ^ 10. The mean comoving distance 
among sources for 10 nearest neighbours is about 150 Mpc. 
Here we compute the mean distance among nearest neigh- 
bours by assuming that all the sources are located at red- 
shift z — 1 and assuming the standard Lambda Cold Dark 
Matter model (|see. for example. Weinber j 12008? ) . Our re- 
sults show that for polarization flux greater than 1 mjy, the 
radio sources show alignment over distance scale of order 
150 Mpc. 

For polarization flux between 0.5 mjy and 1.0 mjy, we 
find a significant signal of alignment for number of nearest 
neighbours greater than 15. The significance is found to be 
greater than 3 sigmas for n v > 30. Hence, in contrast to 
the data with polarization flux greater than 1.0 mjy, here 
we find significant correlations between sources located at 
large distances from one another. We have tested this for 
maximum value of n v = 100, where we find a signficance 
of 3 sigmas. For the entire data sample with polarization 
flux larger than 0.5 mjy we find significant alignment over 
the entire range of n v . At n v = 2 the significance is found 
at 2 sigma level. For all other values of n v the significance 
is comparable or better than 3 sigmas. The distance scale 
corresponding to n v = 100 is roughly 400 Mpc in this case. 

In Fig. 5 we present results for the entire sample of 
12743 sources without imposing any cut. Here we again find 
a statistically significant result for very wide range of values 
of n v . The significance is better than 3 sigmas for the range 
10 ^ n v $C 50. Beyond this range it rises beyond 4 sigmas. In 
the range 150 ^ n v ^ 350 we find a signal more significant 
than 5 sigmas. Hence we see a striking signal of alignment 
over very large distances, of order 500 Mpc, in the entire 
data sample. The signal approaches 2 sigma limit only for 
distances close to a Gpc. In Fig. 6 we show the distribution 
of the statistic Sd for the entire data set for n v — 250. The 
statistic for the observed data is also shown. 

It is also useful to consider the sources which have po- 
larization flux less than or equal to 0.5 mjy. A total of 5291 
sources remain after this cut. The significance of alignment 
as a function of the number of nearest neighbours, n v , is 
shown in Fig. 7. Here we find that sources do not show 
significant alignment over the entire range of values of n v 
shown in Fig. 7. We also probed larger values of n v , up to 
n v = 200. The significance continues to remain less than 2 
sigmas. 

We next consider cuts based on the degree of polariza- 
tion. We expect that the data with low degree of polariza- 
tion would behave the same as the data with low polariza- 
tion flux. However the cut based on degree of polarization 
is not entirely the same as that based on polarization flux. 
In Fig. 8 we show a plot of the total flux as a function of 
the degree of polarization. We notice that the two show a 
strong negative correlation. The Pearson's correlation coeffi- 
cient between the log(Flux) and log(degree of polarization) 
is found to be r = —0.387. This correlation might be in- 
teresting in itself. We also note that the data seems to fall 
into two classes. One shows a clear negative correlation. The 
second with total flux greater than 100 mjy and degree of 
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Figure 3. The statistic So as a function of the number of 
nearest neighbours, n v , for the observed PAs and shuffled PAs 
for polarization flux greater than f mjy (upper graph) The dots 
denote the value corresponding to real data. The points with 
error bars denote the mean value of the random samples along 
with the standard deviation. For larger n v the dots are not 
visible since they merge with the error bars. The lower graph 
shows the corresponding sigma values, i.e. the significance of 
alignment. 
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Figure 4. The statistic So as a function of the number of 
nearest neighbours, n v , for the observed PAs and shuffled PAs 
for polarization flux lying between 0.5 mjy and 1 mjy (upper 
graph). The dots denote the value corresponding to real data. 
The points with error bars denote the mean value of the random 
samples along with the standard deviation. The lower graph 
shows the corresponding sigma values, i.e. the significance of 
alignment. 



polarization lying between 0.01 and 0.1 showing a differ- 
ent behaviour. The degree of polarization is the ratio of the 
polarization flux and the total flux. If the polarization flux 
is uncorrelated with total flux, then the degree of polariza- 
tion would decay like the inverse of the total flux. The data 
shows a weaker decay indicating a positive correlation of the 
polarization flux with total flux. 

We next consider alignment of sources with polariza- 
tion flux less than or equal to 0.01 and those with degree 
of polarization greater than 0.01. For comparison all the 
sources considered in the optical st udy jHutsemekers et alj 
ll998l ; [Hutsemekers and Lamy||2000l ) have degree of polariza- 
tion greater than 0.006. The total number of sources which 
remain after this cut is 3441. We also consider the remain- 
ing sources with degree of polarization greater than 0.01. 
The significance of alignment with these two cuts is shown 
in Fig. 9. Here we find that the data with low degree of 
polarization shows a relatively weak signal for n v =10, 20 
and 30. For larger n v the significance is below 2 sigmas. In 
contrast the data with large degree of polarization shows a 
much stronger signal for n v ^40. Hence we find that sources 
with low polarization flux or low degree of polarization ei- 
ther do not show significant alignment or show a very weak 



signal. The sources with degree of polarization greater than 
0.01 show a strong correlation similar to the one seen in data 
with polarization flux greater than 0.5 mjy. 

We have pointed out above that the data appears to 
fall into two classes, based on the correlation of the degree 
of polarization with the total flux (Fig. 8). We find that 
these two classes show very different alignment statistics. 
The data with total flux greater than 100 mjy and degree of 
polarization lying between 0.01 and 0.1 shows no alignment. 
The remaining data, for which the degree of polarization 
shows a strong (negative) correlation with total flux, shows 
a very strong signal of alignment. There are a total of 10983 
sources in this category. The significance of alignment for 
this set is shown in Fig. 10. The sigma value in this case is 
found to be almost 6 for n v = 150. 

We emphasize that one should be cautious in interpret- 
ing our results for the data sample which includes the sources 
with low polarization flux. The fact that the sample with po- 
larization flux less than 0.5 mjy does not show a significant 
signal of alignment gives us some confidence that the signal 
cannot be attributed to bias. In any case our results provide 
a very strong motivation to make further observations at 
higher precision to further test the signal of alignment. 
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Figure 5. The statistic So (upper graph) and the significance 
(lower graph) as a function of the number of nearest neighbours 
n v for the entire data sample of 12743 sources without imposing 
any cut on the polarization flux. 
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Figure 6. The distribution of statistic Sj), along with a Gaus- 
sian fit, for the complete data set for the number of nearest 
neighbours n v = 250. The statistic for the observed data, 
So = 0.0766 is also shown for comparison. 

5 BIAS 

The fact that low polarizations show very weak alignment 
is very interesting since it strongly disfavors an explanation 
in terms of instrumental bias. If this bias is the source of 
the strong alignment seen in the complete data set, then 
the sources with low polarizations should have shown maxi- 
mal effect. This is because incorrect removal of instrumental 
polarization adds a systematic small polarization to all the 
sources observed in a particular run. For sources with large 
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Figure 7. The significance of alignment for data with polar- 
ization flux less than or equal to 0.5 mjy as a function of the 
number of nearest neighbours n v . 
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Figure 8. The total flux (in mjy) as a function of the degree of 
polarization. 
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Figure 9. The significance of alignment for data with degree of 
polarization less than or equal to 0.01 (plusses) and for degree 
of polarization greater than 0.01 (crosses) as a function of the 
number of nearest neighbours n v . 
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Figure 10. The significance of alignment for data with total 
flux less than 100 mjy or degree of polarization less than 0.01. As 
seen in Fig. 8, this data set shows a strong negative correlation 
between total flux and degree of polarization. 

polarization flux, this will induce a negligible shift in the po- 
larization angle. However the polarization angle of sources 
with low polarization flux will be dominated by this sys- 
tematic effect. Hence we conclude that the effect cannot be 
attributed to instrumental noise and has a physical origin. 

We point out that a bias can generate the alignment 
signal that we observe only if it contributes differently to 
different regions of the sky. If the entire data contains exactly 
the same bias in the polarization angle, then that bias cannot 
yield the alignment signal we observe. This is because, in 
all cases, we compute the statistical significance by shuffling 
the polarization angle among different sources. Hence in this 
case the bias will affect the random samples as much as it 
affects the real data. 

We also perform an explicit simulation to determine 
how the instrumental bias might affect the observed sig- 
nal of alignment. As we have argued above such a bias will 
dominantly affect the sources with low polarization. Hence 
we expect the signal to be much stronger for the data set 
with small polarizations in comparison to the complete set. 
We test this explicitly by a simulation. For this purpose we 
use all the sources in the same positions as the original data 
set. The polarization angles and other observables such as 
intensity, Q and U Stokes parameters are randomized by 
permuting them among the different sources. We consider 
four regions in the sky where the data has instrumental bias. 
These four regions are chosen as (i) RA ^ 6, Dec ^ 30°, 
(ii) < RA sC 6, Dec > 30°, (iii) 10 < RA sC 16, Dec ^ 30°, 
(iv) 10 < RA «S 16, Dec > 30°. Here all the RA values are 
specified in hours. We next induce an instrumental bias in 
this data set by adding a flux of 0.1 mjy to the Q and U 
parameters in regions (i) and (iii) respectively. We also sub- 
tract a flux of 0.1 mjy to the Q and U parameters in regions 
(ii) and (iv) respectively. This mocks the real data which is 
potentially biased due to incorrect removal of instrumental 
polarizations which is expected to be different at different 
epochs of observation. The biased data is then tested for 
alignment. 

After inducing the bias the total number of sources with 
polarization flux less than or equal 0.5 mjy is 5247. The 
number of sources having polarization flux between 0.5 mjy 
and 1 mjy are 3091. The significance of alignment for the 
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Figure 11. The significance of alignment for the simulated data. 
Here we induce an instrumental bias in the data set, as explained 
in text. The upper points (plusses) correspond to the data with 
polarization flux less than or equal to 0.5 mjy. The lower points 
(crosses) correspond to data with polarization flux lying between 
0.5 mjy and 1 mjy. 

simulated data is shown in Fig. 11. We find that the sources 
with polarization flux less than 0.5 mjy show a very strong 
signal of alignment. The maximum sigma value in this case 
is found to be about 9. In contrast the data set with po- 
larization flux lying between 0.5 mjy and 1 mjy shows a 
relatively weak signal. This is clearly in contrast to what we 
observe in the real data. Hence we conclude that the align- 
ment seen in real data cannot be attributed to instrumental 
bias. It is most likely caused by a physical phenomenon, with 
the low polarization sample primarily dominated by noise. 

We have also explicitly verified by simulations that if all 
the data points contain the same intrumental error through- 
out the sky then we do not obtain any significant alignment. 
In this case we find that although the statistic Sd is large 
for the biased data sample, it is also large for the random 
samples generated by shuffling the PAs among the different 
sources. Hence it does not lead to any signal. 

W e next address the issue of bias raised in (|Battve et al.l 
120081 ) . The authors found a significant bias present in the po- 
larization angles in the NVSS survey. It was found that the 
polarization angles have a tendency to be close to multiples 
of 45°. In order to study the existence of such a bias in the 
present data, we plot the distribution of polarization posi- 
tion angles in Fig. 12. Here we show the distribution for the 
complete set, the set with polarization flux smaller than 1 
mjy and the set with polarization flux larger than 1 mjy. 
We find no evidence for such a bias in the present data. 

Finally we study the signal seen in the sample with 
polarization flux greater than 1 mjy for small number of 
nearest neighbours. Here the signal is seen for small number 
of nearest neighbours, n v ^ 14. Since the number of nearest 
neighbours is small it is useful to study the distribution of 
the statistic in more detail to rule out possible large devia- 
tions from Gaussianity. In Fig. 13 we shown the distribution 
of Sd for this set with n v — 10. We find that Gaussian pro- 
vides a good fit to the data. The fit values for 1000 random 
samples are found to be, mean = 0.2836 and standard devi- 
ation^. 00495. If the number of samples are reduced to 500 
we again find a good fit with mean = 0.2834 and standard 
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Figure 12. The distribution of polarization position angles 
(PAs) for the complete data (upper curve) data with polarization 
flux less than 1 mjy (middle curve) and data with polarization 
flux greater than 1 mjy (lower curve). 
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Figure 13. The distribution of statistic So for the data with 
polarization flux greater than 1 mjy and with the number of 
nearest neighbours 10. The random samples are generated by 
shuffling the polarization angles among different sources. The 
Gaussian fit and the statistic for the real data is also shown. 

deviation=0. 00492. Furthermore in this case we can directly 
evaluate the significance of the signal by simulations without 
relying on the Gaussian fit. This determination also agrees 
well with what is obtained by the significance found by us- 
ing the fit. Hence the significance we quote in this case is 
reliable, despite the small number of nearest neighbours. 



6 DISCUSSION 

The possibility that different sources at high redshift 
might show correlations at ver y large distance sca l es wa s 
first indicated by the results of IHutsemekers et al.l (|l99cf ); 
IHutsemekers and Lamvl ()2000h . Here it was found that op- 
tical polarizations from quasars show alignment over very 
large distances. The possibility that such sources might in- 
dicate a simil ar alignment at r adio frequencies was first in- 
vestigated in (|joshi et alj|2007h . They found that the radio 
sources do not show any such alignment. Here we test for 
alignment in radio polarizations usi ng a coordinate inv ariant 
statistics. We confirm the results of ljoshi et al.l (|2007f ) with 



the cuts on the polarization flux that they impose. However 
we do not agree with their conclusions. After imposing the 
same cut as they impose on the polarization flux we do find 
a significant signal of alignment for small number of nearest 
neighbours. The distance scale of alignment here is found 
to be as large as 150 Mpc. If we do not impose any cut on 
the data we find a highly significant signal of alignment for 
a wide range of number of nearest neighbours. The signif- 
icance is better than 5 sigmas for 150 ^ n v ^ 350. This 
corresponds to a distance scale roughly equal to 500 Mpc. 
We also find that the sample with low polarization flux, less 
than 0.5 mjy, shows no alignment. The remaining set with 
polarization flux greater than 0.5 mjy, however, shows a 
strong alignment. Similarly we find that sources with degree 
of polarization less than 0.01 show weak or no alignment. 
In contrast the sources with degree of polarization greater 
than 0.01 show a strong signal of alignment. 

We also find that the degree of polarization shows a 
strong negative correlation with the total flux. The data 
falls into two categories, one of these, which has total flux 
greater than 100 mjy and degree of polarization between 
0.01 and 0.1 shows a random behaviour. The remaining set, 
however, shows a very strong signal of correlation between 
degree of polarization and total flux. We find that this sec- 
ond set also shows a very strong signal of polarization align- 
ment, whereas the first set shows no signal. 

The signal of alignment we find i s similar but may 
not b e identical to that found in l|Hutsemekers et al.l 
1 19981 ; IHutsemekers and Lanwl |2000| ) for optical polariza- 
tions. There the signal of alignment over cosmologically large 
distances was found to be significant for large degree of po- 
larization also. In contrast here we find that if we keep only 
the data with large polarization flux, then alignment is seen 
only over somewhat smaller distances of order 150 Mpc. 
However a very strong signal of alignment is seen over much 
larger distances, of order 500 Mpc, if we do not impose any 
cut on the polarization flux. Furthermore if we relax the po- 
larization cut slightly to keep sources with polarization flux 
greater than 0.5 mjy, we again see a strong signal. Similarly 
sources with large degree of polarization, greater than 1%, 
show a strong signal, where as the remaining set shows only 
a weak signal. 

We also address the issue of whether the signal might 
arise due to bias generated by incorrect removal of instru- 
mental polarization. This bias will dominantly affect the 
sources with low polarizations. Hence if the signal is due 
to this bias then the signal of alignment should be strongest 
for low polarizations. However this is not consistent with ob- 
servations. The low polarization data with polarization flux 
less than or equal to 0.5 mjy shows no signal of alignment. 
The data set with degree of polarization less than or equal 
to 0.01 only shows weak signal for small number of nearest 
neighbours, n v , and no signal for larger n v . Hence it does 
not appear possible to explain the alignment in terms of 
instrumental bias. In any case our results provide a strong 
motivation for testing this signal further with more precise 
data. 

We have not discussed any possible theoretical ex- 
planations for the observed alignment. A potential phys- 
ical effect which could generate such a signal is the 
large scale correlation in the in t ergalactic magnetic field 
(|Agarwal. Kamal and Jainl l201ll ; lAgarwal et afl |2012|) . If 
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the intergalactic magnetic field has a primordial origin 
(|Subramanian and S cshadri 2 0031 ) , then it could show corre- 
lations over very large distance scales. Such large scale cor- 
relations could generate the observed alignment by either 
affecting the electromagnetic radiation during propagation 
through intergalactic medium or by intrinsically aligning the 
sources at high redshift. One possible way in which the back- 
ground magnetic field may affect the radiation propagating 
over cosmological distance s is by mixing with hypotheti- 
cal pseudoscalar pa rticles (|Agarwal. Kamal and Jainll201ll ; 
lAgarwal et alj|2012l ). Pseudoscalar-photon mixing leads to 
very small polarizations at radio frequencies in contrast to 
optical frequencies. Hence the absence of large distance cor- 
relations we found at radio frequencies for large polarization 
flux may be consistent with the predictions of this effect. 
This requires a detailed study which we postpone to future 
research. 
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